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spacecraft surface (thin sheath approximation), nor will it be far 
distant from the spacecraft (orbit*limited approximation), but it will 
be in an intermediate regime such that the— full spacecraft geometry 
must be taken into account. 

TABLE 1. -TYPICAL PLASMA PARAMETERS FOR LEO AND GEO 



GEO 

LEO 

LEO 


(Typical) 

(Low Density) 

(High Density) 

Density (m-3) 

10& 

1010 

1012 

Temperature (eV) 

103 

0.3 

0.1 

Debye Length (m) 

235 

0.04 

0.002 

Electron Thermal Current 
(A/m2) 

10& 

10-* 

10-2 

Ram Ion Current (A/ii|2) 

5x10-10 

10-5 

10-3 

Ion (0*) Mach Number 

0.04 

6.0 

10.0 

Current Collection 

Orbit 

Space -Charge 

Space-Charge 

°CL 

Limited 

Limited 
22 m at 1 kV 

Limited 
3 m at 1 kV 


Currents to a large spacecraft are indicated schematically in fig- 
ure 1. The primary current to the negative surfaces- consists of ram 
ions, while the positive surfaces collect thermal electrons. As the 
ram ion. current is only about one-tenth the electron thermal current, 
the spacecraft will float about 90 percent negative. However, for- high 
bias the ion current will be enhanced by ion.-generated secondary elec- 
trons. An additional source current to negative surfaces is-charge 
blowoff from arc discharges. As discussed elsewhere in this volume 
(refs. 5-7), discharges on solar, arrays have been observed for negative 
biases as low as ^25 0 V. 


y 

PRESHEATH REGION 



Figure 1. A high-voltage spacecraft collects ions 
atnd electrons from a plasma. 
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Positive surfaces Attract plasma electrons, and the resultant 
secondary electrons form a conductive layer over the surface, as dis- 
cussed below. Due to smoothing. of voltage gradients provided by this- 
effect, positively biased array surfaces are much less likely to arc 
than negatively biased surfaces. 

The plasma currents collected by array sutfaces constitute a para- 
sitic current through the array. The effects of this parasitic current 
on. array operation were discussed at the 1978 Spacecraft Charging Tech- 
nology Conference by Domitz and Kolecki (ref. 8). As shown in figure 
2a, solar, cells maintain a nearly constant voltage up to a maximum 
current. When this current is exceeded the voltage drops sharply to 
zero. Clearly, it is advantageous to operate as closely as possible to 
the solar cell's power peak (figure 2b). Figure 3a illustrates that 
plasma current is collected in a distributed fashion over the array. 

The parasitic current at any point of the array is the integral of the 
plasma current density. Figure 3b shows a case in which a relatively 
small parasitic current, when added to the load current, shorts out the 
central portion of the array. 

THE NASCAP/LEO CODE 

NASCAP/LEO is being developed to predict sel£-consistently surface 
potentials, collected, currents, and spatial electric fields for high 
voltage spacecraft in dense plasmas. While much of the code's algo- 
rithmic structure is adapted from the widely used NASCAP code (ref. 9), 
NASCAP/LEO also contains several noteworthy features necessary for 
modeling the intended physical regimes. 

NASCAP/LEC has NASCAP-like object definition routines. Like 
NASCAP, the object is contained in a primary grid (figure 4), which 
may be surrounded by one or more outer grids with successively doubled 
mesh spacing. Enhancements include the ability to put mirror planes 
coincident with one or more of the inner grid boundary planes, and to 
include subdivided regions in the primary grid in order to resolve 
small but important object features. 

Space charge, which is typically ignored in NASCAP, is treated by 
NASCAP/LEO as a local, nonlinear function of potential. Presently the 
function is taken to be 

p - -(e o <P /\ 0 ) 2 [l + m U/a) 372 ]" 1 

A particularly important enhancement in NASCAP/LEO is the ability 
to apply either potential or normal-electric-field boundary conditions 
to each surface cell. The field condition is used for cells whose 
potential is governed by emission of low energy electrons (see below) . 

To calculate plasma currents, NASCAP/LEO defines the sheath sur- 
face as a specified equipotential , and tracks representative particles 
inward. The current represented by each particle is determined taking 
into account ram effects as well as the local plasma density and tem- 
perature. Care is taken not to generate sheath particles in regions 
from which plasma is excluded by high fields or by nearby object sur- 
faces. Thus NASCAP/LEO determines the ion and electron currents to an 
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Figure 3a*. _ Plasma current density, f(V), to a 
solar array (schematic). 
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Figure 3b. Total array current as a function of voltage 

(solid line j schematic), consists of load cur- 
rent, V 0 /Ri , plus parasitic current. For this 
case, a substantial part of the array will be 
shorted out. 


Figure 4. NASCAP/LEO grid structure, showing primary grid (shaded), 
outer grid, and locally subdivided regions. 

object, as well as the distribution of current over the object. 

ELECTRIC FIELD BOUNDARY CONDITIONS 

For electron attracting insulating surfaces dominated by secondary 
electron emission, NASCAP/LEO sets an electric field boundary conditioa 
such that the incident electron current is balanced by a net outflow of 
secondaries to neighboring surfaces. The derivation of electron hop- 
ping conductivity has been presented elsewhere (refs. 10-11), but we 
repeat it here for completeness, then proceed to derive the electric 
field boundary condition. 

Assume a surface has an electron-attracting (positive) surface- 
normal component of electric field, E| , and emits a low energy elec- 
tron current -J e (A/m2) . These electrons travel in parabolic orbits, 
with "flight time", T , given by 

x - 2 [ 2m<e>/e e2j 1/2 

where m is the electron mass and <e> is the (appropriately averaged) 
electron energy. These electrons form a charge density (coul/m2) 
given by -J e t . This- surface charge layer will be accelerated by any 
transverse electric field, E„, to a mean velocity of -e E..T/2m, and thus 
constitutes a surface current, K (A/m) given by 

K * (-J e t2/2)(-e E„/m) = 4 J e E„< e >/E2 = 0|| E„ 

The last relation defines the transverse conductivity, 0 | ( (ohms“*). 

To derive the electric field boundary condition, we further 
assume that the emitted current is proportional to the incident cur- 
rent, Ji n : 



Current balance -requires that 


'in 


En) • 


the previous; derivation assumed uniform fields, it is- no further 
approximation to remove o „ from the divergence, givieg i-urcner 

J in l4( ” Y J in ) < e > /E i] < v * e ii> • 
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ILLUSTRATIVE CALCULATIONS 

i pn illustrate the type of information obtainable from the NASCAP/ 

nf°TflK de P resent . resulCs from three previously published simulations 
of. laboratory experiments. For more details we refer the reader to the 
original publications (refs. 1, 2, 12). co cne 

Early NASCAP/LEO calculations were restricted to potentials about 

Ko^rJJ!i r ^ t f CC i\ eCt:i0n b j’ surfaces of known potential. McCoy and 
ContwIrJ 0 ^ ' 3 ^. re P° rted at the 1978 Spacecraft Charging Technology 
Conference experiments on a biasable metal strip, about 10 meters Iona 

P ! aSti “:' F1 R” 6 5 shows sheath trajectories?^' 
in. which the metal strip was linearly biased. The four cases are 

indicateAi 2 ?h ^A 2 *!* 0 V ’ . and °- 4800 V. The shape of the sheath, as- 

b^lo^lieh? Wo? 1 ? t ? ed . t ^ aject °F Us ’ i& in agreement with observations 

rLJe li 8b 1 5 1 teIeviSTOn - Also, the calculated and measured cur- 

rents-are in good agreement. 

NASCAP/LEO was later improved to predict the spread of hish 
voltages onto insulating surfaces. N. John Stevens (ref. 13) reported 
experiments in which a 3.5 cm diameter biased metal disk (figured) 
was surrounded by either grounded metal or by kapton, and the collected 
current was measured. For the "plain disk" case (figure 7) agreement 

culAa^JhA I he Diak -on-Kapton case (figure 8) is far more diffi- 
cult, as the kapton enhances collected current at high bias voltaee 

?« d f^ , E PreSS6S IC at i?",. v ? lta * e - While the simulation results were 
in fair agreement at high bias, both for collected current and insula- 
tor surface voltage (figure 9) they did not predict the low voltaee 
current suppression. This was because the charging algorithm at the 
time was unable to predict negative charging of tSI kaf?Si sSrfaceS 

A similar experiment was reported by S. Gabriel et al. (ref 151 
who used an emissive probe to measure the electrostatic potential in 
the plasma over a biased "pinhole". The computer model (figure 10) 
made use of the NASCAP/LEO subdivision capability to get good spatial 
resolution in the neighborhood of the pinhole. For negative bUs 

11) agreement was excellent. Fairly good agreement was also 
^f 1 :!^f° r .? 0S ^ 1Ve bias (f ‘« Pre 12). although thl e?p??l.o“ rtoSed 
U£I!i- d i 8 5 inC £ l0 4 between two different plasmas than the calculations 
predicted. Good agreement was dchieved as well for the collected 
currents. This latter calculation showed the utility of the 
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. Current collected by A biased disk surrounded by grounded 
metal. The dashed curve is a prediction by a two-dimensional 
particle-in-cell code (ref. 14). 
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Figure 8. Current, collected by a biased disk surrounded by kapton. 

For high bias, NASCAP/LEO correctly predicts current en- 
hancement (relative to the plain disk). For low bias 
current is suppressed due to accumulation of negative 
charge on the kaptoti. At the time of the calculation the 
NASCAP/LEO charging algorithm did not predict negative 
floating potential for insulators. 






Applied Ground 

Potential 


Figure 9. For high bias, NASCAP/LEO is in good agreement with 
experiment (inset) concerning the spread of high 
voltage onto the insulator. 
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Figure 12. Potentials for 1.27 cm diameter pinhole. 

Solid curves: Calculation. 

Dashed curves: Experiment. 

Upper curves: n e = 2.5 x 107 cm”;?, 0 “ 5.3 ev. — 
Lower curves : n e = 5.8 x 10 cm”-*, 6 “ 4.0 eV. 


subdivision capability, and validated our electric field boundary con- 
dition model for the secondary electron layer. 

FUTURE- DEVELOPMENT 

We plan to develop NASCAP/LEO for use as a reliable design tool. 
Among the planned improvements are prediction of spacecraft floating 
potential, treatment of sheath ionization, improved ram-wake model, 
generalization to full NASCAP geometry, and self-conaistent power sys- 
tem representation. Validation/ improvement cycles will continue at 
each stage, aided by both ground-based and space-based experiments. 
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